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^ Abstract 

q \ In this work we consider the effect of a small-scale helical driving force on fluid with 



a stable temperature gradient with Reynolds number Re < 1. At first glance, this 
system does not appear to have any instability. However, we show that large-scale vortex 
instability appears in fluid despite its stable stratification. In the non-linear mode, this 
instability gets saturated and gives a large number of stationary spiral vortex structures. 
Among these structures there is a stationary helical soliton and a kink of a new type. The 
theory is built on the rigorous asymptotical method of multi-scale development. 

<N . 

^ : 1 Introduction 

\ The importance of the generation processes of large-scale coherent vortex structures in 

hydrodynamics is well known. When these coherent structures appear in small-scale tur- 
bulence they play a key role in transport processes (see for instance PQ). Numerical and 
laboratory experiments [2]- [2] confirm the existence of coherent vortex structures, espe- 
cially for two-dimensional or quasi two-dimensional turbulence [6]-|10j. Notably, they 
are well observed in geophysical hydrodynamics like various cyclones in the planet's at- 
mospheres |12j . The theory of 3D large-scale instabilities is of major interest. An 
example of this theory is the generation of large-scale magnetic field by helical small-scale 
turbulence (vrotv ^ 0) in the MHD (dynamo effect or a— effect) (see for instance |13| . 
[14|). Many works dealt with the dynamo theory generalization for usual hydrodynamics, 
and as a result it was understood that a small-scale turbulence able to generate large-scale 
perturbations cannot be simply homogeneous, isotropic and helical [15], but must have 
additional special properties. So in works [IB] . [T7] it was shown that parity breakdown 
in small-scale turbulence (the external small-scale driving forces) lead to large-scale in- 
stability, the so-called Anisotropic Kinetic Alpha effect (AKA-effect). The injection of 
the helical external force into hydrodynamics systems was considered in many works [18] - 
[21]. In some cases, the existence of large-scale instability was shown (vortex dynamo or 
hydrodynamic a-effect). So, in particular, in work [19] it is shown that the large-scale 
instability exists in convective systems with small-scale helical turbulence. Large-scale 
instability was interpreted as the result of a positive feedback loop between the poloidal 
and toroidal perturbations of the velocity field, which is carried out through the helic- 
ity coefficient. These works, as well as the results of numerical modelling, are in details 



described in review [22J, which is focused essentially on the possible application of these 
results to the issue of tropical cyclone origination. 

In this work, we formulate the problem in a different way. Let us suppose that there 
is a stable temperature stratification in fluid. Let us apply to this fluid with the Reynolds 
number Re <C 1 a small-scale, helical, external force. This force will maintain in the 
fluid small-scale helical fluctuations of velocity field (vrotv 7^ 0). We consider the fluid 
as boundless. At first glance, there are no instabilities at all in this system. However, 
we show in this work that despite stable stratification, a large-scale vortex instability 
appears in the fluid which leads to the generation of large-scale vortex structures. The 
theory of this instability is built rigourously using the method of asymptotical multi-scale 
development similar to what was done in work of Prisch, She and Sulem for the theory of 
AKA-effect [16J . In addition to the linear theory, we also develop and study in detail the 
non-linear theory of this instability saturation. We devote special attention to stationary, 
non-linear, periodical vortex structures which appear as a result of the saturation of the 
found instability. Among these structures there are a spiral vortex soliton and kink of the 
new type. 

Our work is arranged as follows: in Section 2, we set forth the formulation of the prob- 
lem and equations for stable stratification in the Boussinesq approximation; in Section 3, 
we examine the principal scheme of multi-scale development and we give secular equa- 
tions. An overall algebraic scheme of multi-scale development over the Reynolds number 
up to the fifth order is described in Appendix A. In Section 4, we find the velocity field 
of zero approximation and we describe external force properties. Calculation of terms 
of R? order is given in Section 5. The main bulky calculations of Reynolds stress i.e. 
closure of secular equations are presented in Appendix B. In Appendix C we deal with 
some minor questions concerning the closure of temperature equations. In Section 6 we 
consider the overall system of secular equations, and we obtain and investigate equations 
for the large-scale instability. In Section 7, we examine the multi-scale development for 
non-linear cases. The overall algebraic scheme of this development is given in Appendix 
D. In Section 8, we calculate the velocity field of zero approximation for the non- linear 
case. The Reynolds stress calculations for non-linear case are given in Appendix E. In 
Section 9, we discuss the non-linear stage of the instability and its saturation. We study 
the equations of non- linear instability and its stationary solutions. It is shown that due to 
the hamiltonian nature of these equations a large number of stationary vortex structures 
of the spiral type appear. We demonstrate also that there are solutions in the form of the 
spiral soliton and a kink of new type. The obtained results are discussed in conclusions 
in Section 10. 



2 Main equations and formulation of the problem 

Let us consider the equations of motion of non-compressible fluid with a constant tem- 
perature gradient in the Boussinesq approximation: 

dV - - 1 

+ (yV)V = VP + uAV + gpTe z + /„; (1) 

at po 

dT 

— + (V V)T = *AT _ Vz A. (2) 

W = 0,e z = (0, 0, l)-is the single vector in the direction of z axis, /3-is the ther- 
mal expansion coefficient, A = -Q- constant equilibrium gradient of temperature, 
A = Const, A > 0. po = const. VTo = Ae z - the buoyancy force and the external 
force fo,divfo = are taken into account in Euler equation ([I]). Let us write down the 



force /o in the form: /q = /o-^b /~J) where Ao- characteristic scale, to - characteristic 
time, /o- characteristic amplitude of external force. We designate the characteristic ve- 
locity, which is engendered by external force as vo = vo (j^, -^j . When multiplying the 

A 2 

first equation by the parameter we choose dimensionless variables: 

x t -. V P f n P 

Ao t vo Jo poPo 

where 



Ao D ^ , vqv f \l 
v A Ag 1/ 



In the dimensionless variables (t, x, V), motion equations take the form: 

dV f A 2 \ 

— + R(V V)V = -VP + AV + I —5- J ^Te, + F 

^ + i?(FV)T = i-AT - W z (AAo), 
where R = - Reynolds number on the scale Ao, Pr = ~ is Prandtl number. We 

T 
XoA 



introduce the dimensionless temperature T — > t^t and obtain the equations system: 



^ + R(VV)V -AV = -VP + J^-TS, + F , 
1 fdT 1 _\ 



fication we will consider the case Pr = 1. We pass to the new temperature T — > ^, and 



Here Ra = X °^^ - is Rayleigh number on the scale Ao- Further for the purpose of simpli- 
fication we wil 
obtain finally: 



dV 

— + R(VV)V-AV = -VP + RaTe z + F , (3) 
dT \ 

— — AT J = -V z - R(VV)T. (4) 
divV = 0. 

We will consider as small parameter of asymptotical development the Reynolds number 
R = <Cl on the scale Ao- The parameter Ra will be considered neither big nor small, 
without any impact on development scheme ( i.e. outside of the scheme parameters). 

Let us examine the following formulation of the problem. We consider the external 
force as being small and of high frequency. This force drives the small scale velocity 
and temperature fluctuations on equilibrium state background. After averaging, these 
quickly-oscillating fluctuations equal zero. Nevertheless, the non zero terms can occur 
after averaging due to the fact that small non-linear interactions appear in some orders of 
perturbations theory. This means that they are not oscillatory, that is to say of large scale. 
From a formal point of view, these terms are secular, i.e. conditions for the solvability of 
the large scale asymptotic development. So, finding and studying the solvability equations 
i.e. the equations for large scales perturbations, is actually the purpose of this work. Let us 
designate further small scale variables as xo = (xo, to), and large scale ones as X = (X,T). 
The derivative -~ is designated <9j, the derivative is designated d t , and derivatives of 

large scale variables are = V i = df respectively. (No misunderstanding occurs 



between the temperature T and the large scale time T since here time is argument and 
temperature is function). 



3 The multi-scale asymptotical development 

For constructing multi-scale asymptotic development we follow the method which is pro- 
posed in work [16]. First of all, we develop space a nd time derivatives in equations (J3]), 
(|3|) into asymptotical series of the form: 

A = a, + tfv + .... (5) 

^ = 9, + R 4 Sr + -" (6) 
We develop the variables V,T, P like so: 

V(x, t) = v (x Q ) + R(Wi(X) + Vx) + R 2 V 2 + R 3 V 3 + R A V 4 + R 5 V 5 + ■■■ (7) 
Here W\{X)- is the velocity which depends on large scale variables only. 

T(x,t) = T (x ) + R{Qi(X) + Ti) + R 2 T 2 + R 3 T 3 + R 4 T 4 + R 5 T 5 + ■■■ (8) 
Here Q\(X)- the temperature which depends on large scale variables only. 

P(x, t) = ^P-i(X) + Po(x ) + RP 1 + R 2 P 2 + R 3 (P 3 (X) + P 3 )+ (9) 
+R 4 P 4 + R 5 P 5 + ■■■ 

As we will see later, in development of pressure (|9|) it is necessary to have two terms 
which are dependent only on large scales variables P-i(X) i P 3 (X). Let us put now 
the developments ©-Q in the equations system (|3]),([2|) and write down the obtained 
equations up to order i? 5 inclusive. The obtained equations have a rather bulky form and 
are given in Appendix A. In order to simplify the writing of equations we give the algebraic 
structure of development only (vector indices are not written down explicitly, but can be 
easily restored in the necessary places). The conditions of asymptotical development 
solvability ©-([9]) of the equation system ([3]), (jlj) lead to the equations for the secular 
terms (|114p . (|115|) . (|109p . f)105|) and ()104p . Let us write down the full system of secular 
equations: 



3 T W\ - AWf + V p (v p vl + v k Q v p 2 ) + V P (W? W?) = -V k P 3 (X); (10) 



d T @i - A9i = -Vp(t#Zb + v p T 2 ) - V p (WfGi). (11) 
V p Wf = (12) 
W{ = (13) 

VP-i(X) = Ra&i{X% (14) 

The equations (|10p -(|13 p form the main equation system. The equation (|14p is secondary 
and is used to find the field P_i(X), which does not enter in the main system of equations 



(|T0|) - (jT3|) . The line above in equations (|T0j) . (fTTj) means an averaging of corresponding 
terms over quick oscillations. The calculation of these terms (Reynolds stress) is the prin- 
cipal problem. After closing of equations system (|lUp - (|13p it will describe the dynamics of 
large scale perturbations in this scheme. As we can see from equations (jlOp . (jlip . to calcu- 
late Reynolds stresses we have to find fields Vq,v%, Tq, and to average the corresponding 
products over small scale oscillations. In addition, the field v\ is linear depending on W\. 
That is why the Reynolds stresses in the equation (jlOp represent expressions of the form 

Where (i = 1,2) and Pkpi~ the tensor of third rank must be found. The secular 
condition W* = constrain the equation (|10j) very strictly. 
Let us write down the equation (jlOp in components: 

3 T Wf - AW? + V„(W? Wf ) + P xpx V p W? + P xpy V p Wf = -V X P 3 (X) ; (15) 
d T wf - AWf + v p (wf wl ) + /V^Wf + /S^VpW^ = -Vj,P 3 pO; (16) 

PzpxVpW? + PzpyVpW? = -V Z P 3 {X). (17) 

The equation V x Wf + V y Wf = 0, allows us to find from the equations (|15p. ([TUD the 
pressure P 3 (X). But the substitution of this pressure into equation (fT7]) leads to a con- 
tradiction because three equations appear for two variables Wf,W±. There is just one 
possibility to avoid this contradiction, which is to consider the variables Wf, W\ as func- 
tions of variable Z only, i.e. Wf = Wf{Z\W\ = Wf(Z). In this case, the non- linear 
terms in equations (fT5"j) . (fTUj) identically vanish, the equation divW\ = 0, is satisfied 
identically and the equations (fT5|) . (fTBj) take the form: 

3 T Wf - A z Wf + p xzx V z W? + f3 xzy V z Wf = (18) 

3 T W\ - A Z W\ + PyzxVzW? + (3 yzy V z Wf = (19) 
Pz*& g Wf + p ZZ yV z Wf = -V Z P 3 (X). (20) 



Then the velocities are determined by the equations (jT8|h (|19j) . and the pressure is found 
from the equation (|20p . Taking into account the equation (jlip . the temperature 0i must 
also be considered as function of the variable Z only. 

4 Calculations of the zero approximation fields 
(linear theory) 

Let us designate the operator dt — d 2 = Dq. Then, applying this operator to the first 
equation (|1Q2[) . we obtain the equation only for the velocity vo: 

Dy = -Do&Po - Ra{v^l k )r + D F* (21) 

Here l l — vector l l = (0, 0, 1). With help of the equation diV l = 0, we find the pressure Pq: 

PQ = -%w/ lk ^ (22) 



Eliminating the pressure (J22J) from the equation (|2ip . we obtain the equation for vq\ 



DM = -P ip (RaFl k v k ) + D F*, 



(23) 



Here P ip — is the projection operator: 



p%p _ £ip 



QiQp 



As a result, the equation (J23J) can be written down in the form: 

(D$S ik + RaP ip Pl k )v% = D Ft. 
Let us divide this equation by Dq. Then we obtain: 



pip 



(6* + Ra-^Fl"y Q 



k\„.k _ r Q_ 
'o 



Designate operator Lik : 



Pip 

Lik = $ik + Ra-^lplk- 



Then the equation (|26p takes the form: 



pi 

LikV = —, 
Dq 

And the velocity is found using the inverse operator : 

v o ~ L kj Dq ■ 



(24) 



(25) 



(26) 



(27) 



(28) 



(29) 



L ikL kj = 5ij. (30) 
It is easy to make sure by the direct check that the inverse operator L^j has the form: 



L 



RdPkmJ'rrJ'j 
Dq + RdPpqlplq 



(31) 



Consequently the expression for the velocity Vq takes the form: 



RaP^m l m lj 

Dq H~ RdPpqlplq 



Dn 



From the equation 
we can find at once the field To: 

T n = - 



DqTq 



RdP nm l m l n 



Dq ~\~ R&Ppqlplq 



(32) 
(33) 

(34) 



In order to use these formulae we have to specify in explicit form the helical external force 
Tp.The simplest and most natural way is to specify the external force as deterministic. 
(Certainly, it is possible to specify the external force in a statistical way with specifying 
random field correlator, but this leads to more bulky calculations). As it is well known, 
helicity means that forotfo ^ 0. Let us specify the force /o like so: 



/o = fo i cos if2 + j sin ipi + /c(cos 931 + sin <p 2 ) 



(35) 



where 



ipi = k x - u) t, ip 2 = k y - u t, 



(36) 



or 

tp! = klX — OJ t, ip 2 = k 2 x - u) t, , . 

£i = fco(l,0,0); & = *o(0,l,0). 

It is evident that rotfo = k^sfo, where e-is the single pseudo scalar, i.e. helicity is equal 
to : 

forotfo = k ef$ £ 0. (38) 

The formulae f)35j) . f)3T|) permit us to easily make intermediate calculations, but in the 
final formulae we obviously shall take /o,fco,Ct>o as equal to one, since external force is 
dimensionless and depends only on dimensionless space and time arguments. The force 
(|35p is physically simple and can be realized in laboratory experiments and in numerical 
simulation. 

It is easy to write down the force (j35j) in the complex form. It is evident that: 

fo = Aexp(i(px) + A* exp(—iipi) + B exp(iip 2 ) + B* exp(—i(p 2 ), (39) 
where vectors A and B has the form: 

A=^(k- l j),B = ^(J-ik), (40) 
and (fx, ip 2 are given by formulae (|3"7|) . 

5 Calculations of R 2 order terms 

Second approximation equations are the form ()1QT[) : 

D Q v\ = -&P 2 + RaT 2 l l - dk{Wfvi + v k W{) - d k (v k v l + v$v\), (41) 

d t 2 = -v\i k - dkiwfn + v$ei) - dkivfn + v k Tx) (42) 



Let us, as earlier, apply to the equation (I4ip the operator Dq and exclude the field T 2 . 
For v 2 obtain: 

Dlv\ = -D0&P2 - Ra(v%l h )t- 
-RandkiWfn + v k Qx) - d k (v k T + v k Tx)]- 

-Dodk(Wt4 + v k W l x) - D Q d k {v k xvl + v k v\). 

By excluding the pressure P 2 with the help of the equation d{V 2 = 0, we obtain as usual 
the equation with the projection operator (|2"4|) : 



L%u\ + RaP^Fl k v k 2 = 
-P { p {RaFd k (Wx k T + v k Qx) + Ral p d k (v k T + v k Tx)} 
-F* {D d k {W k vl + v k Wf) + D d k (v k v p + v k vl)\ . 



We divide this equation by D 2 , and see that it can be written down in the form: 



pip 
Pi 

rk a ,P 



{RaPd k (W?T + ugeO + RaPd k (vfT + w*T a )+ 



(43) 



Where L^- is the former operator (|27j) . Taking into account the expression (I3ip for the 
inverse operator L7. , we obtain: 



D%+RaP Xq l x l q 



{RalPd m {WpT + v^@ 1 )+ 



(44) 



+RaPd m {v?T Q + v^T x ) + D d m (W^ + «) + £> a m ««g + «)}• 



Since t>i, Ti^o, To-are functions of only fast variables xo,then the second and the fourth 
terms in Reynolds stresses Vt>ot>2 do not give contribution in formula (|44p as well, we can 
omit them. Thereby, we will be interested in the expression: 



4 



RQ j PkrJ j rJ j j 

D 2 + RaP Xq l x l q 



pjp 

-2 (RaFWrdmTo + D Q W?d m v p ) 

i 



(45) 



o 



Let us write down the expression (|45p in the form: 



„,k _ jx/-inrpmk Ti/mrpmk 



Where the tensors Tffi,Tj%F have the form: 



(46) 



ink 



(1) 



?kj 



RaPk n l n lj 
D 2 + RaP Xq l x l q 



pjp 

-—y-Ral p d m To, 
D o 



(47) 



rpmk 
1 (2) 



R(lPkn^n^j 

D 2 + RaP Xq l x l q 



pop 

7T'" 



-d m Vn 



(48) 



The calculations of Reynolds stresses are rather bulky, which is why we perform them in 
Appendix B. 



6 Large scale instability 



In Appendix B, we calculate Reynolds stresses for the equation (jlOp : in other words, 
we obtain closed secular equations for large scale perturbations. We write down these 
equations in the explicit form: 



d T W x + aV z W y = V 2 z W x , 



(49) 



d T W v - aV z W x 



VtWy 



(50) 



_ 4-2Ra 

a = —eRa- oxo 51 

(4 + Ra 2 ) 2 v ' 

Here e designate a single pseudo scalar because expressions VzW y , —VzW x are compo- 
nents of rotW. The equations ([4"9~]). (f50|) differ from equation of AKA -effect [16] by the 
coefficient a only. They obviously contain an instability which generates large scale vortex 
structures. Choosing the velocities W x ,W y in the form: 



W x = AexpfaT) s'mkz, 



(52) 



Wy = B exp(7T) cos kz, (53) 

we obtain the instability increment 7 = ±ak z — kg, i.e. max 7 = with the k = 
The formulae ( l52jh (f53|) describe spiral vortex structure (circularly polarized plane wave) 
with the amplitude which increases exponentially with time. These waves are sometimes 
called Beltrami runaways since for them there is no usual hydrodynamical interaction 
WVW = 0. 

If the external force has zero helicity, then the a-term vanishes in accordance with 
the general theorem of the Reynolds stress tensor [15]. Helicity is taken into account in 
the external force structure itself. If the temperature gradient vanishes, then it is evident 
that the a-term also vanishes. Besides the equations (|49p . (|50p there is an equation to 
find the pressure: 

C (W x + W y ) = -P 3 (X) + Const, (54) 

where C 

_ 8 - \2Ra - 2Ra 2 - 3Ra 3 - Ra A 
C ~ (4 + i?a2)3 ' (55) 

and also the equation for the large scale temperature perturbation found in Appendix C: 

dr^ - AGi = -2 (4 ^°~ a 2 )2 Vz (W x + W y ) (56) 
and the equation to find the pressure P-i(X) : 

V Z P^(X) = Ra0 1 (X) (57) 

It is easy to solve all these equations with the known fields W x , W y . 

It should be noted that the field Wf is equal to zero in the main approximation 
only. W z 7^ in the approximation of higher orders. In compliance with this in the 
approximation of higher orders appear also derivatives Jjj?, ^7. This means, that in fact 
one has to consider that 

d d d 

dz^dx'W (58) 

The conditions (I58p means that the horizontal scales of the instable vortices L x , L y ^>L z 
(vertical scale). L z 3>Ao- 

Since the increasing fields W x ,W y belong to Beltrami type, the general non linearity 
(VFV)iy cannot limit their increase, because, as it has been already stated, it is identically 
equal to zero. In order to describe the process of this instability saturation it is necessary 
to develop a non linear theory of a-effect. This theory is described in the following 
sections. 



7 Multi scale development for the non linear case 

In the non linear case the large scale field W(X) is no small any longer that is why the 
asymptotical development (HI)-© has to be modified. Let us search the solution for the 
equations ([3]), (j4j) in the following form: 

V(x, t) = 4#_iP0 + v (x ) + RV X + R 2 V 2 + R 3 V 3 + ■■■ (59) 
R 

T(x, t) = 4 T -iP0 + T o(x ) + RTi + R 2 T 2 + R 3 T 3 + ■■■ (60) 



P(x,t) = ^P^(X)+^P^ 2 (X)+^ 1 (X)+P (x )+R(P 1 +P 1 (X))+R 2 P 2 +R 3 P 3 +- ■ ■ 

(61) 

For the derivatives we use the previous development ©,([6]). Substituting these expressions 
into the initial equations (|3|) , (|3J) and gathering together the terms of the same order R, up 
to degree R 3 inclusive, we obtain the equations of multi scale asymptotical development. 
The algebraic scheme of this development is given in Appendix D. The result of this 
scheme is the separation of secular terms from fast oscillating ones. Let us give secular 
terms in the explicit form: 



d T W l - AW 1 + V k v k vl = -ViPi; (62) 



d T T - AT + V k (v k T ) = 0. (63) 

In these equations we do not write the law index (—1). Besides this, there are secular 
equations: 

VfW* = 0, W z = 0, (64) 

V k {W k W % ) = -ViP-i, (65) 

V k (W k T) = 0. (66) 
The equations ([64"|) - (|66j) are satisfied in the previous geometry: 

W = {W x (z),W y (z),0),and,P- 1 = Const. (67) 
There is also an equation to find the pressure P_3 : 

V Z P_ 3 = RaTl z . (68) 

It is clear that the essential equations for finding the non linear alpha-effect is the equation 
([62]) , In order to obtain these equations in the closed form we need to calculate the 



Reynolds stresses Vfc [vqVqJ . Below, we will deal with the solution of this problem. 

First, we have to calculate the Vq fields of zero approximation. From the asymptotical 
development in zero order we have the equations: 

d t vl - d 2 vi + W k d k vl = -diP + RaToP + i& (69) 
d t T -d 2 T + W k d k T = -v k l k . (70) 

8 Calculation of zero approximation fields in the 
non linear case 

Let us introduce the operator Do : 

D = d t - d 2 + W k d k . (71) 
Using the operator Dq, we write down the equations ([69]) and ([70]) in the form: 



D v l = -d t P + RaT l l + F \ 



D T Q 



-v k l k 



Eliminating the temperature from the equation (|72|) we obtain: 

D 2 v = -D diP - Ra(v k l k )r + D Q F*. 

dA = 0. 

Eliminating the pressure from the equation (|74p we obtain: 

kjki 



or: 



D£v l = -P ip (Rav*l k P) + D Fq 



(D 2 5 ik + P ip Ral k F)v% = D F*. 



Dividing this equation by Dq, we can write it in the form: 

T ,v k - ^0 
Do 



(72) 
(73) 

(74) 
(75) 

(76) 
(77) 

(78) 



Where is the prior operator (|27p in which Do is replaced by Do- In much the same 
way the inverse operator also coincides with ([3T]) when substitung Do by Do- As a result 
it is easy to find the fields Vq and Tq: 



RaP] im l m lj 
^^-^pg^p^ci 



Ra,P nm l m l n 

^^Ppq^p^q 



Do 



(79) 



(80) 



The external force Fo has the prior form (|39p . The effect of the operator Do on proper 
function exp(iu>i + ikx) has obviously the form: 

Do exp(iu;t + ikx) = Dq(uj, k) exp(icji + ikx), where Dq{u, k) is: 



From this it is evident that: 



D (cj,k) =i{u + kW) + k 2 . 

D (ui, -fa) = i(u - faW) + kf, 

Dq(u, -ki) = D (-u>, fa), 
D (u, -k 2 ) = i(oj - k 2 W) + k%, 
D* (u,-k 2 ) = D (-u:,k 2 ). 



(81) 

(82) 

(83) 
(84) 
(85) 



From the formulae (|79p and (|39p follows that the field Vq is composed of four terms:t>Q 

v 01 + v 02 + v 03 + v 04 where 



v 02 — (^Ol)*) ^04 



v i — e 



Skj 



RaPk m l m lj 



Dl{-ujo,fa) + RaPll 



A> 



D (-u)o,fa) 



(86) 



J 03 



RciPk m l m I'j 



D 2 (-uj ,k 2 ) + RaPll 



D (-u ,h) 



As it was said earlier, in scalar operators Dq one can take ujq = 1, fci 
(0, 1,0). Then taking into account that Pll = 1, we obtain: 

D (u ,-k 1 ) = l + i(l-W 1 ) = D 1 , 



(87) 

(l,0,0),fc 2 = 



D (-u> ,k 1 )=D* 1 



D (cj ,-k 2 ) = l + i(l-W 2 ) = D 2 , 
D (-u} , k 2 ) = D 2 .. 



(89) 



Taking into consideration these formulae we can write down the velocities Vq in the form: 



"01 



4s 



Skj 



RaP] tm l m lj 
Df + Ra 

RaP] tm l m lj 
Df + ifo 



P.J 



(90) 
(91) 



9 Non linear instability and large scale vortex 
structures 

The calculations of the Reynolds stresses for the non linear case are performed in the 
Appendix E. Let us write down in the explicit form the equations for the non linear 
instability: 



d T Wi - VlWi 



(92) 



Ra(l - W 2 ) 



[1 + (1 - W 2 ) 2 } [(W 2 (2 - W 2 ) + Ra) 2 + 4(1 - W 2 ) 2 } ' 



-V, 



d T w 2 - v 2 z w 2 = -v.rg 

Ra(l - Wi) 



(93) 



[1 + (1 - Wi) 2 ] [(Wi(2 - W x ) + Ra) 2 + 4(1 - Wx) 2 } ' 

Here we introduced the following notations: W\ = W X ,W 2 = W y . It is easy to verify 
that with small values of the variables W\ , W 2 the equations (l92j) , (|93j) are reduced to the 
equations (|49p -(|5ip and describe the linear stage of instability. It is clear that with the 
increasing of W\, W 2 the non linear terms decrease and the instability gets saturated. As 
a result of the development and stabilization of instability the non linear vortex helical 
structures appear. The study of the form of these stationary structures is of interest. For 
that purpose we take 8tW\ = 8tW 2 = in the equations ([92]) . (f93|h Integrating these 
equation over z, we obtain: 



RaP 



dX _ 

dz ~ (1 + P 2 ) [AP 2 + (1 - P 2 + Ra) 2 } 



dP 

dz 



RaX 



(l + X 2 )UX 2 + (l-X 2 + Ray 



C 2 



(94) 
(95) 
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Figure 1: Phase picture of the dynamical system with Ra — 2, C\ — — 0. The bold line 
shows the phase trajectory which comes out of the point (1,1) and after the "time" Z = L 
comes back to the same point. This trajectory presents the stationary solution of the boundary 
problem with the rigid boundaries in the layer whose thickness is L = z. 



Here new variables are introduced X = 1 — W±, P = 1 — W2, and C\, C2- are integration 
constants. The system of equations (|94|) and (|95p can be write down in the hamiltonian 
form: 



dX _dH dP _ dH 
Ih ~ ~dP , ~dz~ ~ ~dX' 
Here the variable z plays the role of time and the hamiltonian H has the form: 

H = U(P) + U(X) + dP + C 2 X + C 3 . 
Where U (x) has the form: 



(96) 



(97) 



U(x) 



1 



■In 



2\2 



/Ra 



4(4 + Ra) ARa + (x 2 + 1 - Ra) 2 4(4 + Ra) 



arctan 



1 + x 2 - Ra 
2y/Ra~ 



(98) 



The function H (J9T|) . (i98l) is obviously the first integral of the equations system (i94"|) . (1951) 
and can be found by the direct integration of this system. With C\ = 0, C2 = the 
function U (x) is limited above and below as well. That is why the hamiltonian section 
by the constant H = Ho, gives closed periodical trajectories on the phase plane (X,P), 
which correspond to the helical vortex structures in the real space. Examples of phase 
pictures for Ra = 2 and Ra = 3 are represented in fig. 1 and fig. 2. With C\ = 0, C2 = 
on the phase plane there is only one elliptical point. Closed trajectories correspond to the 
periodical non linear vortex structures. Thick closed lines correspond to the non linear 
structures which are also the solutions of the boundary problem with the rigid boundary 



W x = 0, W 2 = 0, z = 0, z = L, 



where L is the period over z of phase trajectory, which gets out with z = 0, of W\ = 
0, W 2 = and gets back to the same point with z = L. The space structures of periodical 
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Figure 2: Phase picture of the dynamical system with Ra — 3, C\ — 0, C2 = 0. The bold line 
shows the trajectory which corresponds to the stationary solution of the boundary problem 
with rigid boundaries with z = and z = L. 



solutions is presented in fig.3-fig.5. If one of the constants, for instance C\ 7^ 0, then one 
hyperbolic point appears on phase picture. For instance, phase pictures with C\ = 0.1 
are presented in fig. 6. The example of periodical vortex structure which corrersponds to 
the closed trajectory on pase plane with Ra = 2 is given in fig. 7. The solution which 
corresponds to the separatrix in fig. 6 is of particular interest. This solution describes the 
solitary spiral turn of the velocity field around the axis z (soliton) see fig 8. Moving away 
from soliton the velocity field becomes constant. This kind of solitons were not known 
earlier. The interesting particularity of this soliton is the fact that it is also the solution 
of the boundary problem with free boundaries. For this boundary problem [23j : 

dWi _ dW2 _ Q 

dz dz 

on the fluid boundary. In addition to that, boundaries must be on a big distance from 
the soliton, much bigger than the soliton's characteristic dimensions. In cases when there 
are two constants C% 7^ 0, C2 / two hyperbolic and two elliptical points appear on 
phase picture. The example of this phase picture with C\ = 0.1, C2 = 0.1 is shown in 
fig.9. As earlier, the periodic vortex structures correspond to closed trajectories around 
elliptical points. Localized solutions (solitons) correspond to a separatrix in fig.9. Since 
the separatrix connects two different hyperbolic points, the soliton now has two different 
limiting values, with z — > ±00 , fig. 10. This soliton is called a kink. Thereby, spiral kinks 
correspond to the separatrix in fig.9. These kinks are also solutions of the boundary 
problem with free boundaries. In conclusion, it should be remembered that the system of 
the equations (l92j) . (1931) is closed. The velocity field Wi,W% determines the pressure Pi, 
according to the formulae: 

~r\ — 1^) -t- 1( 2 ), 

where , Tg are given by the formulae (|235[) . ()236[) . Besides, the velocity field Wi, W2 
gives the contribution to the equation for temperature ()63p . Closure of this equation is 
made in much the same way as the closure for velocity. Nevertheless, this equation is 



Figure 3: Helical vortex structure with Ra = 3, G\ = 0, C 2 = 0. 



secondary and here we do not give the result of this closure. 

10 Conclusions and discussion of results 

In this work, it is shown that in fluid with stable stratification, a large scale instability 
appears under the action of small scale helical force. The result of instability is the gen- 
eration of vortex structures of the Beltrami type. The vortices have the characteristic 
vertical dimension L z ^>\q and the horizontal dimension is much bigger than the verti- 
cal one. Since the vertical component of the velocity W z is equal to zero in the main 
approximation and the stratification is stable, then the found instability has no relation 
with convection. The structure of the equation which describes the instability in linear 
approximation is the same as the equation of a-effect or more precisely as the equation of 
AKA-effect. As a result, instability generates plane spiral waves with circular polarisation 
(Beltrami runaway). With an increase in amplitude, the instability and its stabilization 
are described by a non linear theory. Stationary equations appear to be hamiltonian, 
which is why they are a rich set of periodical spiral vortex structures. Notwithstanding 
the fact that attention in this work was essentialy paid to a boundary free problem, it 
should be noted that some periodical solutions turn out to be solutions of the boundary 
problem with rigid boundaries. We would like to pay special attention to stationary soli- 
ton and kink, which correspond to the separatrix on the phase plane. This is the solitons 
of the new type. In real space it describes one spiral turn of the velocity vector field 
around the axis z. Soliton and kink are also the solutions of a boundary problem with 
free boundaries. 

Let us return to the formulation of the problem. The external helical force /o is given 
in the explicit form in order to make calculations more transparent. Strictly speaking, 
its explicit form is not very important for the existence itself of a-effect. It is necessary 
that rotfo~f only. The external force could be chosen statistically by specifying the 
correlator : 

fif m = A(t, r)5 im + B(r, r)rir m + G(r, r) (99) 

It is fundamental that the last term G(r, r) (helicity) in this correlator is not equal to 
zero, otherwise a-effect is absent. Nevertheless the statistical method is more bulky since 
it requires the specification of the functions A,B,G and calculations of rather compli- 
cated integrals. If we specify the external force dynamically then the averaging over fast 



Figure 4: Helical vortex structure with Ra =1.8, C\ = 0, C 2 = 0. 



oscillations is performed easily. 

The question is of interest about the origin of instability on qualitative level. For this, 
we revert to the expression for the Reynolds stresses: 

V z (v z v§ + v$vl). (100) 

As far as the direction z is particular, the averages (v^v^ + VqV^) which are non equal 
to zero must be proportional to l z . From the property of the external force ^.l^Uj = 0, 
follows, that with k = 1, the velocity W x i fully left of Reynolds stress. Another velocity 
W y enters into Reynolds stress with the coefficient iW y l z B^Bl. Since from the properties 
of external force it follows that B3BI = — |, we obtain the factor (jW y lz). In a similar 
way, with k = 2, because of property of the external force -B2-B3 = 0, the velocity W y 
is fully missed out in the second Reynolds equation. The velocity W x enters on the 
second Reynolds equation with the coefficient: iW x l z A\A^. Due to the property of the 
external force ^4^3 = f> we obtain the factor (—\W x l z ) . (We may not write down the 
common factor). It is clear that as a result we obtain the components of vector product 
[W x /] = iW y l z — jW x l z . Or if we take into account V z ,we obtain the components 
rotW. These components rotW provide a positive feedback loop between the velocity 
components like in the usual a-effect which leads to instability. 

Appendix A. Asymptotical development scheme 

In the i? _1 order there is only one term : 

ap_i(x) = o. (101) 

With this the equation (jlOip is satisfied automatically. The equations of the order R° has 
the form: 

d t v - d 2 v = -dP + RaT l z + F , 

d t T - d 2 T = -ug, (102) 
dv Q = 0. 

From the equation f|102[) it follows immediately that the functions vq,Tq, Pq are oscillating 
due to the oscillating character of the external force Fq . Approximation equations R 1 have 
the form: 



Figure 5: Example of a vortex structure with Ra = 1.6, G\ = 0, C 2 = 0. 



d t vx - 8 2 Vl + dv v = -8P 1 - VP-x{X) + Ra(Qi(X) + T x )l z - 

dtTi - d 2 T x = -W{ - vf - dv T ; (103) 

dvi = 

The equations (I103P contain already oscillating terms as well as the non-oscillating ones. 
Oscillating terms after averaging give zero and only the non-oscillating ones remain. That 
is why the solvability condition of this system is the independent vanishing of the non- 
oscillating (secular) terms as well as the oscillating ones. For the system f|103[) the condi- 
tion of solvability in the approximation R 1 gives equations: 

VP-i(X) = RaQ x {X)l z (104) 



Wf = (105) 
The oscillating part in the approximation R 1 is described by the equations system: 

d t vi - d 2 vi + dv v = -dP\ + RaT x l z \ 

d t T x - d 2 n = -v{ - dv T ; (106) 
dvi = 

i.e.vi = v\(xo),Ti = Ti(xo),P\ = P\{xq). In the approximation R 2 we have the following 
equations: 

d t v 2 -d 2 v 2 = -dP 2 - d(WiVQ + voWi) - d(viVQ + VQVi) + RaT 2 l z ; 
d t T 2 -d 2 T 2 = - v *-d(WiTo + v Qi)-d(v 1 T + v T 1 ); (107) 
dv 2 = 0. 

It is easy to evidence that after averaging, all the terms in equations (|107j) give zero. 
Thereby secular terms do not appear in the order R 2 and the fields v 2 ,T 2 ,P 2 remain 
oscillating. However, they now depend on large scale variables X, i.e. v 2 = v 2 (xq,X),T 2 = 
T 2 (xq,X), P 2 = P 2 (xo,X). The approximation R 3 gives equations: 

dm - d 2 v 3 = -8P 3 - 8(Wi + vxXWx + Vl )- 

-d(v 2 v + v v 2 ) + RaT 3 l z ; , , 

8 t T 3 -d 2 T 3 = - v ^-d(W 1 +v 1 )(e 1 + T 1 )-d(v 2 To + v T2); [ ' 

dv 3 + VWi = 0. 
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Figure 6: Phase picture of a dynamical system with Ra = 2, Ci = 0.1, C 2 = 0. 

From the equation (|108p one can see that an averaging of the first two equations gives 
only zero terms, i.e. does not give any secular terms. But the third equation averaging 
gives a new secular term: 

VWi = (109) 

Thereby the fields v 3 ,T 3 ,P 3 remain oscillating, but depend on X, i.e. v 3 = v 3 (xq, X),T 3 = 
T 3 (xq, X), P 3 = P 3 (xq, X)- Now we pass to the equations of the i? 4 approximation. They 
have the form: 

d t v A - 2dVv 2 - d 2 v A + d[v 3 v + v v 3 + v 2 {Wi + vi) + (W% + v 1 )v 2 }+ (110) 

8 t T 4 - 2<9VT 2 - 8 2 T A = -vl-8 [v 3 T + v T 3 + v 2 (@i + Ti)+ 
+(Wi + Vl )T 2 ] - V[(Wi + vi)T + wo(0i + Ti)]; 

dv 4 + Vv 2 = 0. (Ill) 

It is easy to see that these equations have no secular terms at all. 

Finally, let us consider equations of the approximation^ 5 . Their form is rather bulky: 

d t v 5 + d T Wi - d 2 v 5 - 2dVv 3 - AWi + d[viV + v Q v 4 + v 3 (Wi + v%)+ (112) 
+(W 1 +v 1 )v 3 +v 2 v 2 }+V[v 2 v +v v 2 +(W l +vi)(W 1 +v 1 )} = -dP 5 -VP 3 (X)-V P 3 +RaT 5 l z . 
8 t T 5 + a T 9i - 8 2 T 5 - 28VT 3 - AGi = -v z 5 - d[v A T + v T A + v 3 (&i + 7i)+ (113) 



+(Wi + vx)T 3 + v 2 T 2 ] - V(v 2 T + v T 2 ) - V[(Wi + vi)(9i + H)]. 
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Figure 7: Helical vortex structure with Ra = 2, Ci = 0.1, C 2 = 0. This structure corresponds 
to the closed trajectory around the elliptical point in flg£J 



dvs + V^3 = 

The equations of the fifth order f)112j) . (|113p give the main system of secular equations as 
the solvability conditions of the fifth approximation: 



d T W 1 - Alfi + V(v 2 v + v v 2 ) + V(WWi) = -VP 3 W; 



(114) 



d T @i - A6i = -V(v 2 T + v T 2 ) - V(Wi9i 



(115) 



Appendix B. Calculations of the Reynolds stresses. 

Since the external force is composed of four terms (|39p . then the expression for Vq is 
composed of four terms as well: ±Vq, 2 Vq, 3 Vq,± Vq. With this : 



2V = (lV ) ,4 V = { 3 V ) 



(116) 



iv Q = expttpi 



3 v Q = expup 2 



RaP pX (h)lxh 



D%(-u , h) + RaP mq {ki)l m l q 



A,, 



RaP p \(k 2 )l\l^ 



D%(-u , k 2 ) + RaP mq (k 2 )l m l q 



D (-uj ,ki)' 
By 

D (-ujo,k 2 )' 



(117) 



(118) 



In order to simplify equation writing, we will write down the convolution P mq l m l q m 
the formP(A:)//. In the similar way for the temperature field To there are four terms: 
iT),2 Tq,3 Tq,4 Tq. At the same time : 



2 Tq — (iT )*,4T) — (3T0)" 



(119) 



1^0 



- exp upi 



RaP(h)ll 



D%(-u ,ki) + RaP(k 1 )U 



(120) 
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Figure 8: Helical soliton which corresponds to the separatrix in fig. 6 with Ra — 2, C\ — 0.1, 
C 2 = 0.1. 



3-^0 



exp iipi 



RaP{k 2 )U 



D%(-u ,k 2 ) + RaP(k 2 )ll 



Dl{-u ,k 2 ) 



(121) 



Let us consider the tensor T^-y 

rpmk rpmk rpmk r pmk. 
11 (2) > 2i (2) ' 3J (2) ' 4J (2) • 



It is also composed of four terms: 



rpmk / rpmk\* rpmk / rpmk\* 

2 J (2) - U J (2) ) >4 J( 2 ) - U i (2) ) ■ 



122) 



The equation for the tensor follows directly from the formula ()48p : 



rpmk 
1 T (2) 



exp up i 



n '■ ] 



1 D%(-uoM)+RaP(kl)ll 



P„(k!)ik m 



D^(-uj ,k 1 )+RaP(ki)ll 



(123) 



A„ 



First of all we transform the expression, which contains projection operators in the formula 
(fT23|) . 



RaP kn (ki 



D%(-u; ,ki) + RaP(k 1 )ll 



P 



jp 



RaP p \(ki)lxlf, 



Dl{-UQ,ki) + RaP{ki)ll 



RaP kn (k)l 

n<"i 



D%(u ,k) -RaP(k)ll 



P 



Raid 



Dl(uj ,k) - RaP(k)ll 



(124) 



Pkn 



5n\ 



RaP\j(k)l n lj 



D%(uj ,k)+RaP{k)ll 



Ral\L 



D$(u} ,k) + RaP(k)ll 



Here we use the property of the projection operator. 

PjpPpX = PjX 

With help of the formulae ()124p it is possible to write down: 
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Figure 9: Phase picture of the dynamical system with Ra = 2, C\ — 0.1, C 2 = 0.1. One 
see the appearance of two hyperbolical and two elliptical points. 



Tmk "'"-i 
- (2) =expz^i— g 



ik?P kn (h) 



D^(-L0 ,kl 

8\n ~ 



RaP\j{ki)l n lj 



DK-ojqM) + RaP{kx)ll 
Ralxlfj, 



(125) 



DK-uq,^) + RaP(h)ll 



ink 



. ik?P kn (k 2 ) 
(2) =expz( / 9 2 — T — 



RaP X j(k 2 )l n lj 



1\n 



D%{-uj ,k 2 ) + RaP{k 2 )ll 
Ralxl^ 



(126) 



Dl(-uo,k 2 ) + RaP{k 2 )U 



Bp. 



From the definition ()47|) follows expressions for four terms of the tensor: 

(127) 



rpmk. rpmk rpmk rpmk rpmk . 

' 2 J (l) '3 J (l) ' 4 " t (l) • 



■ink 
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mk \ * 



I rpmk\* rpmk / rpmk\ 
U-t(l) J ,4-t(l) - (3-t(i) J 



rpmk 
lT (l) 



ik^RaV 1 ^ 



^-expz^i-^- 

RaP(k x )ll 
D%(-u ,ki) + RaP(h)ll 



Pkn(h)>< 



(128) 



.4 



(1) 



ik^RaW 



-u , k 2 



-Pkn(h)x 



(129) 



RaP(k 2 )ll 



D%(-Lu ,k 2 )+RaP(k 2 )ll 



B„ 



As a matter of fact, all these expressions are considerably simplified. Actually, the opera- 
tor Dq = dt — d 2 acts on its own function exp(iu;oi + ikx) and gives Dq ex.p(iuot + ikx) = 
Dq{ujq, ko) exp(iuot + ikx), where Do(uo, ko) = iooo + /eg. In dimensionless variables this 
means that : 
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Figure 10: Helical kink which corresponds to the separatrix in figEJ 




D (uj , k ) = 1 + i; D (-u , k ) = 1 - i = Dq(uj , k ). 
The expression Pknhln = Pll = l 2 — ^jS- = 1 5 since ki,k 2 -i-l- As a result: 

Do (w , ko) + i?aP// = 2i + Ra; D 2 (-oj , k ) + RaPll = (2i + fla)* 
And all tensors are simplified considerably: 
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(130) 
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The reason of the further simplification of these expressions is : 



Ppu(k)l u Ip. 



Taking into account (|138|) . we obtain: 
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We will do calculations of Reynolds stresses in several stages. To begin with, we consider 
the term v^T^. This average value is composed of four terms in which the oscillation 
phase is cancelled: 



v P T$ = [ lV P] [ 2 T- fc ] + [ 2 <] [ X T$] + [ 3V P] [ 4 T™ k ] + [ 4 <] [ 3 T™ k ] (145) 

The second term in the (|145p is conjugated with the first one, and the fourth with the 
third one. Now we substitute in the (I145P correspondent expressions for tensors and we 
obtain: 
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RaFA z 
(2i + Ra)* 



(146) 
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RaHuBl 



(1 + *) 
>w the 

consideration the (|146p . we obtain: 



Ral p B z irp . _ , fcR , _ 2Ral k B* z 

(2t + Ra)* l[ k 2 y (2i + Ra) ^ (2i + Ra) 2 ' 



Let us find now the components of Reynolds stresses v^v k = —VqW^T^. Taking in 



~^p~k o iW x 

V ° V2 = Re (TTi) 
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(2i + RaY 



(147) 



A* — h k A* 
A k K l A x 



2Ral k A* z Ra%A* z 
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■\B P 



Ral p B z k 2Ral k B* 



(2i + .Ra)* 



Ra 2 l k B* 



(2i + Ra) (2i + Ra) 2 > 



The full contribution in the tensor of the Reynolds stresses fgwf + v^v^ from the tensor 
Tjffi is obtained using the symmetrization of this equation over the indices p, k. As a 
result we obtain: 



P k i £• P r> *"as 

VqV% + VqU^ = -Re 
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(148) 



»W y i?aP5 2 ir „, lfc „, 2Ral k Bl 

a+o [ 



Ra 2 l k B* 



(2t + Rci)* 11 k 2 y (2i + i?a) (2* + i?a) 2j 



+ 



(1+0 



L4 A 



Ral k A z 
j2iTRaj 

Ral k B 



-][At-k{Al 



2Ral p A* z 



+ 



Ra 2 LA 



2i + Ra) ' (2i + i?a) 2 ' 
2RaLBt Ra 2 LBl 



+ 



.(1 + (2i + #a)* J[ p 2 « (2i + i?a) ^ (2i + Ra) 2 ^ 

Let us put p = z, k = 1, in the equation (|148p . i.e. we find the component x in the 
Reynolds stress. It is easy to see that: 



V V 2 + V V 2 



ReV 2 



+ReV z 



iW y Ra 
(1 + ~ (2i + Ra)* 

2Ra Ra 2 



+ 



(1 + (2i + Ra) (2i + Ra) 2 



]B z Bl \ + 
}B* z B l 



Taking into account the B* Z B\ 



4' 



we obtain: 



V l V 2 + V Q V 2 



(149) 



(150) 



ReV z W v ; 

y \4{l + i) 



Ra 



2Ra 



Ra 2 



(2» + i?a)* (2* + Ra) (2i + #a) 2 
After calculating the real part in the (|150p we obtain: 



Where : 



v o v 2 + v o v 2 = aiV z W y , 

„ (-4fla + i?a 2 + 12) 
-ita- 



(151) 



(152) 



4(i?a 2 + 4) 2 

Putting p = z, k = 2, in the (11481) . we find the corresponding component in the Reynolds 
stress. It is easy to see that : 



(1 + 

^ a a A * 2Ra % Ra 2 % 

(U + RaY^ 2 ~ j2JTRa) A2A3 + J2iTRa) 2 2 3 



(153) 



+ 



2Ra 



+ 



Ra< 



(2« + Ra) (2i + i?a) ; 
As far as = 0,^3^2 = f> we obtain from the (|153[) : 



B2BI 



vlv\ + ugvf = -aiV^. (154) 

Now we need to calculate the contribution in the Reynolds stresses from the tensor 
Iff. As it was done previously 



„.Vrp m k 
V 1 {1) 



[ 2 T ( Tf ] + [2<] [iT ( -f] + [ 3 <] UT$] + [ 4 <] [ 3 T ( ^] (155) 

With this the second term is conjugated with the first one and the fourth is conjugated 
with the third one. The simple calculation of the first term gives: 

2 



Ra 



1 - 



Ra 



Ap 

We calculate similarly the third term: 



2(1 + i) 3 (2i + Ra) 

RalpAs 



-ik^)l k x 



(2i + Ra)*_ 



A* 



1 - 



Ra 



2(1 + i) 3 (2i + Ra) 

RalpBz 



-ik?)l k x 



B, 



(2i + Ra)* 



Bl 



Now it is easy to find the contribution in the v^v^- 



Re 



-Re 



Ra 



(1 + i) 3 

Ra 
~(l + i) 3 



1 - 



Ra 



V b 
V V 2 = 

2 



(2i + Ra) 



Ra 



(2i + Ra) 



(iW x )l l 



(iWy)l h 



A, 



B, 



Ral„A 



(2i + Ra)* 

RalpB^ 
(2i + Ra)* 



Bt 



After symmetrizing this tensor over the indices p, k, we obtain: 



V k 1 k P 
V V 2 + V V 2 



(156) 



(157) 



(158) 



(159) 



Re< - 



Ra 



+Re 



+Re 



+Re 



H + if 

Ra 
{l + if 

Ra 
(1 + i) 3 

Ra 

"(TmF 



Ra 



(2i + Ra) 



Ra 



(2i + Ra) 

Ra 
(2i + Ra) 

Ra 
(2i + Ra) 



(iW x )l h 



(iW y )l h 



(iW x )P 



{iWy)F 



A,, 



B n 



At 



Bi 



Ral„A 



p^3 



(2i + Ra)* 

RalpB^ 
(2i + Ra)* 

RahA 3 
(2i + Ra)* 

RahB 3 
{2i + Ra)*_ 



A%\ + 



Bt > + 



A\ 



Bt 



Putting p = z,k = 1, in the formula (|159p . we obtain the tensor x-component of the 
Reynolds stress: 



V V 2 + V V 2 



(160) 



„ iRa 
Re< -- xrr 1 



Ra 



(1 + i) 3 V {2i + Ra) 
Since AiA% = 0, BiB^ = |, we get: 

Ra 



Vqv\ + v\v\ = Re ■ 



. 4(1 + f) 3 

After calculating the real part, we obtain: 



V z [W x AiA* 3 + W y 5i^] 



(2i + Ra)J VzWy 



Where 



V V 2 + V V 2 = OSVjgWj, 

(4 - Ra 2 - ARa) 
° 2 ~ ~ Ra A{A + Ra 2 ) 2 
Putting in the formula (|159p p = z,k = 2 we easily get: 



«o v 2 + V U 2 



Re I — — I 1 



Ra 



(1 + i) 3 V (2i + Ra) 



V z [W X A 2 A* 3 + W y B 2 B* 3 ] 



Since A 2 A* 3 = -\,B 2 B% = 0,then 



v o v 2 + v o v l = -a2V 2 W x 



and designating : 



a = a\ + a 2 = —Ra 



A-2Ra 



1611 



(162) 



(163) 



(164) 



'(4 + i?a 2 ) 2 ' 

we obtain the final expressions for the tensor components of the Reynolds stresses: 



v o v 2 + v o v 2 = otV z W y , 



(165) 



(166) 



(167) 



v%v\ + vlv z 2 



-aV z W x . (168) 

The terms (|167p and (| 168j) are fundamental. Nevertheless, strictly speaking one must 
calculate other, less important, terms. Let us consider the component p = z, k = 3 of the 

pmk 

-(2) 



Reynolds stress. Putting in (|148p p = z, k = 3, we get the contribution of the tensor T™ fc- 



VqV 2 + l$vf 



(169) 



i?eV 2 



+ J ReV, 



2tW g 
T+i 

2iW y 
T+i 



Ra 



(2i + Ra)* 

Ra 

(2i + Ra)* 



2Ra 



+ 



Ra' 



(2i + Ra) ' (2i + Ra) 2 
2Ra Ra 2 



(2i + Ra) (2i + RaY< 



A 3 A* 3 J. + 



B 3 BI 



As far as ^3^3 = -B3-B3 = \, we obtain: 



where 



ugvf + wjjvf = C 2 V Z (W x + W y ) , 
-Ra 3 + 2Ra 2 - 4Ra + 8 



(4 + Ra 2 Y 



(170) 



(171) 



Now we can find the contribution in VqV^ + v^v^ of the tensor T™^. Putting in the formula 
(I159j) p = z, k = 3, we obtain: 



V0V2 + v o v l = Re ^z < 7; 



Ra 



(2i + -Ra)* 



(W X + Wy) 



v z vl + v 3 v* 2 = dV z (W x + W y ) , 



where 



Ci = -Ra- 



2 + Ra 



(4 + Ra 2 ) 2 

Designating the common coefficient C = C\ + C 2 , 

8 - 12Ra - 2Ra 2 - 3Ra 3 - Ra 4 



C 



(4 + Ra 2 Y 



we obtain: 



v z v 3 + vfr% = CV Z (W x + W y ) 



(172) 
(173) 

(174) 

(175) 
(176) 



11 Appendix C. Closure of the temperature 
equation. 



In order to close the temperature equation we have to calculate the term: 



V z [vf 2) T {0) + vf Q) T {2) 



(177) 



It follows from the formula (|42f) , that indeed there is the contribution only of the following 
terms in Tr 2 )' 



■(2) 



(178) 



'0 v ' ^0 

From the (|34j) it is easy to find the component iT( ),2 ^~(o)>3 ^(o)j4 ^~(o)> i^Tq) = 2^(0) j 



(0) -4 i (0) 



2 T (0) 



4 T (0) 



exp(-i(/9 2 ) 



(1 + i) 2 
1 

JT+W 



Ra 



(2i + Ra) 

Ra 
(2i + Raj 



A* 
^3> 



Do 



1 



T (2) = WfT^ + WfTf'T - —Wpd p T { 



W D 



(o) 



(179) 
(180) 
(181) 



First of all we find the component V 2 f^T( ) : 



V*«f 2) T (0 ) = 2ReV z { [^f 2) ] [ 2 T (0) ] + 



3«(2) 



[4%] } • 



(182) 



Since the u?^ is composed of two terms (|46p . then at the beginning we find the contribution 
-W" 9 [iT^][2T( )]. As it was done in previous calculations we obtain: 



W q [iT$}[ 2 T {0) ] 



-i)RaW x 



(l-i)4(l +i )2 
Further we find the contribution: 



Ra 



(2i + Ra)* 



Ra 



-W%T^][ 2 T {0) ] 
As a result we get the term: 



Ra 



2ReV 2 



2i + Ra)* 
[2%)] 



(2* + Ra)* 



Ra 



A Z A* 



(2i + Ra) 



A Z A* Z . 



(2) 



(183) 



(184) 



(185) 



2ReV z W x {--) 



Ra 



1 



1 



Ra 



(2i + i?a)" 



4 y -if 

Let us come now to the calculations of the term 2ReV z 

2ReV z <, :i r 



Ra 



(2i + Ra) 



A Z A* Z . 



{ [3^2)] U T (o)] } : 



(2) 



(o)J 



= 2i?eV, {-W g [ 3 2j*][ 4 r (0) ] - W ff [ 3 2g][4r ( o)]} 
Calculations similar to the previous ones give : 



-i)i?aW y 



(l_i)4( 1 + i )i 



-^[3^][4T (0 )] = ^ 



^[3Tg][ 4 r {0) ] 

_ Ra 
~ (2i + Ra)* 

Ra 



Ra 



1 



(2* + Ra)* 



(2i + i?o) 
Ra 



B Z B* Z , 



1 



(2« + fta)* 



As a result we obtain the term V 2 f^T( ): 



V ^(2) T (0) 



(186) 



(187) 



(188) 



(189) 



2ReV z W x ( 



Ra 



{I -if 



Ra 



(2t + ifa)* 



(2i + i2a) 



|AJ 2 + 



+2i2eV z W v (- T ) 



(1-i) 5 



Ra 



(2i + Ra)* 



Ra 



(2» + i?a)* 



I5 2 | 2 . 



Let us come now to the calculations of the term u? %T(2). 



V*wf 0) T (2) = 2Re ■ 



1^(0) 



[ 2 T (2) ] + 



3» 



Let us write down the auxiliary expressions: 

1 



[2 T (2)] : 
[4 T (2)] 



(1 + 
1 



(190) 

(191) 
(192) 



Taking into account these formulae and expressions for [2??^], [2^^], [4^'^], [4^ J, it is 
not difficult to get an expression for the V z i^T( 2 ) : 



(0) J (2) 



-2ReV, 



iW x 



2(1 + i) 2 



Ra 



{l + if 



Ra 



Ra 



(2i + Ra) 
2 



+ 1 



Ra 



(2i + Ra) 



+ 1 



(2i + Ra) 



2ReV, 



Ra 



{l + if 



2(1 + i) 2 

Ra \ ( Ra 
+ 1 



(193) 



(2i + Ra) J V ( 2i + -^a) 
2 



+ 1 



1 



Ra 



\B, 



(2i + Ra) 

The sum of the expressions (|193p and (|189p is the Reynolds stress (|177p . After simple 
algebraical transformations and finding of the real part we obtain the final expression: 



Vz [vf 2) T {0) + vf Q) T {2) I l J ^_l_ Vz {Wx + ^ 



(4 + Ra 2 ) 2 



(194) 



12 Appendix D. Scheme of asymptotical develop- 
ment for the non linear case 

Let us present the algebraical structure of the asymptotical development of the equations 
©! @ f° r the non linear theory (we will not write indices because they can be restored 
evidently at any moment). In the order P -3 there is only the equation: 



&P_3 = 0,^P_3=P- 3 (X). 
In the order R~ 2 we have the equation : 

ap_ 2 = o,=^p_ 2 = p_ 2 (x). 

In the order P _1 we get a system of equations: 



(195) 
(196) 



d t W-i - d 2 W-i = -(0P_i + VP_ 3 ) + RaT_ x l z - dW-{W-x, 



(197) 



d t T_x - d 2 T_! = -dW^xT-x - Wi x , (198) 
dW-x = 0. 

The system of equations (|197p . f|198j) gives secular terms: 

- VP_3 + RaT-xlz = 0. (199) 

W z _ x = 0. (200) 
In zero order R° we have the following system of equations: 

d t v - d 2 v + d{W^xv + v W„x) = (201) 

= -(dP + VP_ 2 ) + RaT l z + F, 
8 t T - 8 2 T + d(W^xT + v T^x) = ~v z . (202) 
dv = 0. 

These equations give one secular equation: 

VP_ 2 = 0, =>• P-2 = Const. (203) 
Consider the equations of the first approximation R: 

d t vx - d 2 vx + d(W-xVx + v x W-x + v v ) = (204) 

= -V(W-iW-i) - (3Pi + VP-i) + RaTxh. 

d t Tx - d 2 Tx + d(W-xTx + Wl r_i + « 7b) + V(W-xT-x) = -vf. (205) 

dVx + VVF_i = 0. (206) 
From this system of equations follow the secular equations: 

VW-x = 0, (207) 
V(W-xW-x) = -VP_i, (208) 

V(W_iT_i) = 0. (209) 

The secular equations f|20T[) - (|2Q9j) . are obviously satisfied in the previous velocity field 
geometry: 

W= (W x (Z),W y (Z),0);T„x = T_i(Z);VP_i = 0,=> P-i = Const. (210) 
In the second order R 2 , we obtain equations : 

d t v 2 - d 2 v 2 - 2dVv + d{W-xv 2 + v 2 W^x + v vx + vxv ) = (211) 



= -V(W-xv + v W-x) ~ {dP 2 + VP ) + RaT 2 l z , 



8 t T 2 - d 2 T 2 - 2dVT + d(^_iT 2 + v 2 T„ 1 + v Ti + Vl T ) = 



(212) 



= -V{W- 1 T + v T^ 1 )-v 2 . 

dv 2 + Vv = 0. (213) 

It is easy to see that in the order R 2 there is no secular terms. 

Let us come now to the most important order i? 3 . In this order we obtain equations: 

d t v 3 + drW-i - (d 2 v 3 + 2<9Vt>i + AW_i) + V(W_itq, + uiW_i + v v )+ (214) 

+d(W- 1 v 3 + W3WL-1 + u U2 + v 2 v + = ~(dP 3 + VPi) + PaT 3 Z z . 
d t T 3 + 9 T T_! - (<9 2 T 3 + 23VT! + AT_i) + V(W_iTi + mT_i + v T )+ (215) 

+9(iy_iT 3 + « 3 T_i + u T 2 + v 2 T + «iTi) = -v z 3 . 

dv 3 + Vwi = 0. 
Prom this we get the main secular equation: 

d T W-! - AW-! + V (vpo) = -VPi, (216) 
S(rT_i - AT_i + V(w To) = 0. (217) 

13 Appendix E. Reynolds stress in non linear 

case 

In order to calculate the Reynolds stresses we have first of all to calculate the expression: 



v*vh = 2Re(y* 1 v&+v* 3 v&). C2\s) 
Taking into account the formula (|90j) . we obtain: 

1 

\Di\ 



« + « = T$ = 7^ (A k A* + A* k Ai)+ (219) 



RaA* z ( l k Ai + UA k \ RaA z ( l k A* + 1 { A 



|^| 2 V D 2 + Ra J \D X \ 2 V Df + Ra 

2 Ra 2 l k U \A Z \ 2 
|£>i| 2 \D 2 + Ra\ 2 ' 
Similarly taking into account the formula (|9ip . we obtain: 

1 

\D 2 



+ 



<4<3 + <4Xs = T$ = TT^iBkB* + BIB,)- (220) 



RaB* i l k B{ + liB k \ RaB z f l k B* + hB k ( 



D 2 \ 2 \ Dl + Ra J \D 2 \ 2 V Df + Ra 



the components of the tensor T?, 



| 2 Ra 2 l k k \B Z 
\D 2 \ 2 \Dl + Ra\ 

(1) and T (2) 



It is clear that the components Tf\ and T?X are of interest. To begin with we consider 



(I)" 

1 

\Di\ 



T (i) = 7^(^*1 + (221) 



SinceA^ = A%A X = 0. 



Ra ( A\A X A 3 A* X \ _ 

\D\ + Ra + Df + Ra ) ~ ' 

T (i) = 77^(^2 + ^4 2 )- (222) 
1^1 1 

Ra ( A\A 2 , A 3 A* 2 \ 



+ 



|Di| 2 \Dl + Ra Df + RaJ 
The first bracket in the ()222|) is equal to zero, which is why: 



, 32 _ i Ra {Dl-Df) 

(i) ■ 

Now consider the component T?2 : 



jM ' ^ -1 -1 ) (223) 

W 4 \ Dl \ 2 \Df + Ra\ 2 

■(2) ' 



T (2) = TTT^ 3 ^* + 5|S 2 )- (224) 

1-^21 



Ra 



\D 2 \ 2 \D 2 + Ra Df + RaJ 
As far as B^B 2 = B^B^ = we consider the component T^y 

1 

\D 2 



T ( 3 2 ) = TW^^ B * + B * Bl) - (225) 



i?a / B*B X B 3 Bl \ 



\D 2 \ 2 \Dl + Ra Df + RaJ 
The first bracket in the formula (l225p is equal to zero, then: 



, 31 = . Ra {Df-Dl ) 
{2) ~ ^\D 2 \ 2 \Dl + Ra\ 



m = -7 ,^,2, J „ (226) 
Taking into account : 

(D 2 - Df) = - Wi), (Df - D 2 ) = -4?(1 - W 2 ), (227) 

\D X \ 2 = 1 + (1 - Wi); \D 2 \ 2 = 1 + (1 - W 2 ), (228) 

\D\ + .Ra) 2 = (Wi(2 - Wi) + Ra) 2 + 4(1 - Wi) 2 , (229) 

\Dj + Ra\ 2 = (W 2 (2 - W 2 ) + Ra) 2 + 4(1 - W 2 ) 2 . (230) 
The components T?j 2 ,, take the form: 



T 32 = Rajl-W^ 

« [1 + (1 - W^W^ - Wi) + Ra) 2 + 4(1 - Wtf] ' 

T 3i = Mi - w 2 ) 

(2) " [l + (l-W 2 ) 2 ][(W 2 (2-W 2 ) + Ra) 2 + A(1-W 2 ) 2 }' 
Now calculate the components and Tj^y It is easy to see that: 



t-i33 _ 



rpS3 

1 (2) 



2|A 3 



2 IB, 



1 - Ra 



1 - i?a 



Ra 2 



{Dl + Df) 

\Dj + Ra\ 2 \Dj + Ra\ 

(D 2 + Df) Ra 2 



Or in the explicit form: 



^33 



\D 2 + Ra\ \D 2 + Ra\ 



1 



1 -Ra 



Ra 



(i) 2 [1 + (1-W 1 ) 2 ] 

2[1 - (1 - Wi) 2 ] + Ra 
[Wi(2 - Wi) + i?a] 2 + 4(1 - Wi) 2 

rp33 — v 

V~ 2[1 + (1-W 2 ) 2 ] X 

2[1 - (1 - W 2 ) 2 } + Ra 
[W 2 (2 - W 2 ) + Ra] 2 + 4(1 - W 2 ) 2 



(231) 
(232) 

(233) 
(234) 

(235) 



(236) 
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